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Abstract 

The method of unitary clothing transformations put forward in relativistic quantum field theory 
(QFT) by Greenberg and Schweber and developed by Shirokov is applied to construct a new family 
of interactions in the meson-two-nucleon system. Along with a brief exposition of its basic elements 
we show a specific transition from the initial "bare" one-meson and one-nucleon operators and 
states to their physical "clothed" counterparts. We emphasize that the clothing transformations in 
question do not alter the original total Hamiltonian, but provides a conceptually more transparent 
representation of the same Hamiltonian in terms of a new set of operators for particles with physical 
properties and their relativistic interactions. The Hermitian and energy- independent interaction 
operators for the processes vrA^ — > ttN, NN — > A^A^ and A^A^ <-> ttNN are derived starting from 
the Yukawa-type couplings between fermions (nucleons and antinucleons) and bosons (vr— , 77—, 
p—, cv— mesons, etc.). These types of interaction have a distinctive off-energy-shell structure 
which is naturally generated by the unitary transformation that removes from the Hamiltonian the 
(three-leg) vrNN vertex coupling. 
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I. INTRODUCTION 



To view the light hadronic systems as few-body systems involves several problematic 
aspects, and many of these aspects still need to be theoretically explored in detail. Not 
all cases are so fortunate that one can accurately reduce the underlying field-theoretical 
problem to a non-relativistic quantum-mechanical Hamiltonian formulation in a (two-body, 
three-body, etc.) potential scheme. One fortunate situation occurs typically in low-energy 
few-nucleon physics ^|J, where the considered degrees of freedom are the nucleons, and the 
non-relativistic interactions are generated via meson-exchange processes, or via alternative, 
effective field-theoretical derivations. But, even in this well-defined research area, there have 
been attempts to brake the paradigm of the purely nucleonic Hamiltonian description in the 
potential scheme Q, and to relate the origin of well-known unsolved problems, such as the 
Ay puzzle, to dynamical effects of the pion 0. 

As the energy increases, the light hadronic systems become much more problematic to be 
described. The opening of the pionic threshold forces to deal directly with the pion degrees of 
freedom. This complicates the few-body equations since one has to deal simultaneously with 
phenomena of production and absorption of particles Q . The derivation of non-perturbative 
dynamical equations for meson-few-nucleon processes is still waiting for a consistent theory 
that solves the general problem in a feasible and calculable way. 

Specific topics, such as the study of NN — > A^A^vr inelasticities, both in free space and 
in more complex nuclear environments, have received much attention recently. Many phe- 
nomenological works have attempted the description of the experimental data, including 
polarization data, and extensive reviews on the past and current status of research on this 
subject are available 0, S, i, Q H. In general, such models perform some nonrelativistic 
reduction of the interaction Hamiltonian between mesons and nucleons, and construct the 
basic production mechanisms using perturbative schemes. Typically, these terms consist 
of the ttNN one-body term (direct term plus recoil) and some more complex two-nucleon 
processes such as pion rescattering, intermediate A excitations, and 2N processes involving 
heavy- meson exchange currents. As yet, there is not a clear-cut understanding of which 
mechanisms contribute the most and, especially in the isoscalar channel at threshold (the 
channel of the pp— ^-ppvr" process) there has been quite a debate on which are the basic 
mechanisms that provide the strength to the threshold cross section. This strength can be 
generated either by short-range heavy- meson exchanges or by pion rescatterings in the 
a channel Presently, the analysis of the A^A^ — > NNti inelasticities did not succeed to 
resolve this ambiguity, and there is hope that the study of more complex reactions, such 
as pd 7T° He^, can help to solve this ambiguity. Phenomenological studies for this reac- 
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tion have been limited to calculations involving only pion rescattering in the a-channel 
these calculations achieved, very recently, a good reproduction of data for both angular 
distributions for the excitation function and for the analyzing power [12]. For this 3N reac- 
tion, a thorough study involving production mechanisms originated also by 2N heavy-meson 
exchanges is needed, but unfortunately these studies have still to come. 

Besides the problem of constructing the 2N heavy-meson exchange mechanisms in a 
consistent manner (we will come back to this problem in the final part of this paper), there 
remain, in addition, several difficulties that need to be explored. One is related to the 
need to use relativistic equations for the nucleon dynamics in the energy range of meson 
production, since this energy range is higher than the typical low-energy nuclear regime 
where nonrelativistic quantum mechanics works so well. Then, the first step is to construct 
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consistent relativistic (pseudo) potentials that fulfill the Poincare algebra. 

Another difficulty is connected with the treatment of virtual effects generated by the 
ttNN vertex in the production process. To directly use the one-body vertex and similar 
operators in the construction of one-body-type production mechanisms may create serious 
difficulties because in the nNN vertex the three particles cannot stay simultaneously on- 
mass-shell. However, the use of quantum-mechanical equations for the nucleon dynamics 
forces the nucleons on their mass shell, and the pion, being the physical observed particle 
in the production process, has also to stay on its mass shell. Clearly, all this creates am- 
biguities in the possible treatment of one-body operators, and hence there are problems of 
internal inconsistencies in the construction and calculation of sets of one-body, two-body, 
etc, operators. This represents a not negligible problem in pion production at threshold: 
although the magnitude of the one-body term is generally small, the interference effects of 
the one-body term with respect to the more dominant two-body amplitudes have important 
consequences in the cross sections. 

A similar difficulty also occurs when the elementary fermion-meson vertex interaction is 
used for the construction of meson decays operators in relativistic constituent quark mod- 
els. It has been shown recently that the use of the spectator-model ansatz (built upon 
the elementary quark-meson vertex) in Relativistic Quantum Mechanics (RQM) is not fully 
constrained in covariant calculations of strong (meson) decays of baryon resonances. This 
aspect has been pointed out in Ref. [3| specifically for the construction of spectator-model 
operators in the so-called Point-Form formulation of RQM. However, the presence of ambi- 
guities in the general construction of Current Operators has been observed in all forms of 
RQM, including front form and instant form ^ . 

The approach illustrated in this paper represents a field-theoretic way to overcome the 
above mentioned difficulties. Among the varieties of field-theoretical approaches we concen- 
trate on the Unitary- Transformation (UT) methods for Hamiltonian-based models. UT's do 
not change the expressions for the observables, because the S-matrix remains the same under 
such transformations. However, the unitary transformations allow to change the represen- 
tation of the Hamiltonian so that the interaction operators refers to particles with physical 
(observable) properties. It is important that the interactions extracted this way are Her- 
mitian, energy-independent and include the off-shell structures naturally. These attractive 
features of the unitary transformation approaches allow to derive pion-production interac- 
tion operators that are relativistic and involve only physical particles (and hence, one-body 
operators are eliminated from the representation and all the virtual processes or effects are 
embedded in the interaction operators) ; we also extend the calculation to include 2N heavy- 
meson exchange mechanisms for the production amplitude, given the importance that these 
mechanisms have in discussing the phenomenology of 3N forces and pion production. 

There have been quite different implementations of UT techniques (see Ref. [1^ and refs. 
therein, for a review). In such approaches one starts with the original field Hamiltonian H, 
to arrive at a new representation of H which is unitarily equivalent and where the particles 
involved interact through interaction operators that are energy independent. Such energy 
independence simplifies the solution of the eigenvalue problem for the initial Hamiltonian. 
To avoid the confusion that already exists in the literature, we do not employ here the 
inflated term "effective" to qualify our interactions, since they have no direct connections 
(that we are aware of) with the formalisms of effective Lagrangians. 

Amongst the different UT formalisms, we mention first Okubo's idea to block- 
diagonalize H via a UT with respect to the decomposition of the full Fock space IZp oi 
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hadronic states into the nucleoli (no- meson) subspace and its complement in TZp. The 
Okubo approach was further developed by Glockle and Miiller their relativistic the- 

ory of interacting particles. These authors were the first to show, al least to our knowledge, 
how with the help of one single UT the non-commuting Hermitian operators (Hamiltonian 
and Lorentz boosts) can be reduced to a block form. For a simple model of "spinless nu- 
cleons" exchanging scalar mesons, it was then possible to construct the new Hermitian and 
energy independent contributions to the Poincare group generators. 

In Ref. [l8| the same method was employed for deriving the nucleon-nucleon and nucleon- 
antinucleon interactions starting from a field Hamiltonian with the exchange of tt— , p— , uj— 
and a— mesons. Then, in the framework of the Hartree approximation, these interactions 
have been introduced to describe the saturation properties of nuclear matter. That approx- 
imation, we note in passing, gives rise to an original recipe (see Eqs. (23)- (25) in [18]) for 
calculating the nucleon mass shift in the three dimensional approach. 

Fuda and co-workers [l^, jl^], [1^, l^l have also used the Okubo-Glockle-Miiller ap- 
proach to construct one-boson-exchange models in light-front dynamics. They prefer to work 
with the mass-square operator rather than with H because the former commutes with 
all of the Poincare generators. Such studies demonstrated that the two-particle interactions 
obtained in second-order perturbation theory have certain similarities but do not coincide 
exactly with the relevant Born terms obtained using the Feynman-diagram techniques (see 
also Ref. [l3)- We shall address this question in more detail in Sec. 3. 

Successful applications of the UT method have been obtained also in the theory of nu- 



clear forces and nuclear electromagnetic processes (see Refs. j23[, [2J] with reference to 
earlier studies of the 50's, e.g., Ref. [1^). Such developments by the Osaka group jl^l il- 
lustrate another useful application of the UT formalism: the elimination from H of certain 
virtual processes whose S-matrix elements, by definition, have to vanish on the energy shell. 
Remarkably, the transformed Hamiltonian H' = WHU contains only interactions for real 
processes. Even more remarkably, it was possible to derive the one-pion-exchange (OPE) and 
two-pion-exchange (TPE) nucleon-nucleon potentials, the pion-nucleon potential (see also 
[2^), the NA —> NN transition potential, the interactions for pion absorption/production, 
and so on, from one and the same initial Hamiltonian H without any additional ad hoc as- 
sumption (of physical or mathematical nature) in each case. The corresponding Hermitian 
energy-independent operators are defined in the complete Tip, i.e., not only in a subspace 
C Tip as was the case for Okubo approach. 

In the present work, we construct the meson-baryon interaction operators relying upon 
the notion of " clothed" particles. Such notion was introduced many years ago by Greenberg 
and Schweber ji^l (see also Chapter XII in the monograph j^) in their aim to include the so- 
called cloud or persistent effects (the terminology refers back to Van Hove [2^, [1^). In this 
connection, we apply the method of unitary clothing transformations elaborated by and 
[3^ and developed recently in Refs. ji^, js^], js^, LL5i]- In this way, a large amount of virtual 
processes induced by the meson absorption/emission, the NN -pair annihilation/production 
and other cloud effects can be accumulated in the creation (destruction) operators for the 
clothed (physical) mesons and nucleons. Such a bootstrap reflects the most significant 
distinction between the concepts of clothed and bare particles. 

The approach used in this article differs from the aforementioned ones at least in two 
aspects. First, the clothing procedure does not aim, a priori, to define a UT that blockdi- 
agonalizes H, as is done in the Okubo approach. Besides, it might be not always possible 
to implement such a task in the infinite-dimensional Hilbert space. Instead, the aim of 
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the multistep clothing procedure is to express the original Hamiltonian H in terms of the 
new clothed-particle operators in a form which is different than that given for the ini- 
tial bare-particle ones. Such a transition from the bare-particle representation (BPR) to 
the clothed-particle representation (CPR) introduces a new sparse structure of the original 
Hamiltonian H. Second, in the framework of the CPR the mass and vertex renormalization 
problem 3|, js^ is considered in a natural way, in parallel with the construction of the 



interactions, while the renormalization problem in Ref. j2J| has been disregarded. Other 
distinctions are discussed in Sects. 2 and 3. 

In this context, we would like to draw attention to the recent work by Stefanovich on the 
problem of divergences in QFT. As was shown in Ref. [s^, one can introduce a similarity 
transformation WHU to cancel infinite counterterms directly in H and therefore also in 
the ^-matrix operator. As an application, the method has been illustrated for the case of 
an Hamiltonian system in quantum electrodynamics. Although Ref. [s^l does not use the 
clothing ideas, there are many points of contact between that approach and the CPR method 
and, in our opinion, further developments in this area could have promising perspectives. 

The central goal of this paper is the derivation of interaction operators in the CPR for 
the processes ttN — >■ ttN, NN — > NN and A^A^ ^ irNN. For the field-theoretical treatment 
of mesons and nucleons we assume the instant form after Dirac of relativistic dynamics, in 
which the generators of space translations and rotations are the same as in the free theory. 
We remind that in this case the three boost operators must involve the interaction parts to 
meet, together with the total Hamiltonian H, the well-known commutation relations of the 
Poincare-Lie algebra. At the beginning we prefer to consider if as a function of creation 
and destruction operators in the Fock space rather than the more customary form that 
follows from a given field Lagrangian. Our option follows a general statement that H may 
be expressed as a sum of products of creation and destruction operators (see Chapter HI of 
Ref. |38|). This Hamilton formulation of RQFT simplifies the introduction of the CPR that 
we consider in this work. 

The outline of this paper is as follows. The aspects of the method of unitary clothing 
transformations, which are necessary for constructing a new family of interactions operators 
in the meson-two-nucleon system are illustrated in Sec. 2. We introduce an auxiliary UT's 
that convert the primary bare bosons and fermions into some intermediate- level particles 
with physical masses. The approach develops along the chain: bare particles with bare 
masses — *• bare particles with physical masses physical (observable) particles. This pro- 
cedure is useful for drawing some parallels between the clothing approach in QFT and the 
method of canonical transformations (in particular, the Bogoliubov ones) in the theory of 
superfluidity and superconductivity. Analytical expressions for the relativistic interactions 
(quasipotentials) in the CPR are shown in Sec. 3 in case of the primary Yukawa coupling 
between pions and nucleons. Our calculations are performed using purely algebraic means 
within an iterative technique proposed in that Section. The features that distinguish the 
obtained momentum-space quasipotentials from their on-energy- shell counterparts derived 
in the second-order Dyson perturbation theory with Feynman rules are emphasized. Ex- 
plicit formulae for the quasipotentials generated by the additional heavier-meson exchange 
mechanisms are given in Sec. 4. 
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II. FORMALISM 



A. Introductory definitions 

Our departure point is the Hamiltonian 

H = H{a)=Ho{a)+Vo{a), (1) 

where the unperturbed Hamihonian HQ{a) and the interaction term Vo{a) depend on the 
creation and destruction operators of the "bare" particles with unphysical masses and cou- 

o 

phng constants. Here, a denotes the set of all these operators. For example, in case of a 
spinor (fermion) field ip and a neutral pseudoscalar meson field one has to introduce the 

o o o 

operators a (k), b (p, d {p,r) and their adjoint counterparts, respectively, for mesons, 
nucleons and antinucleons. These operators enter in the expansions 

(x) = (27r)"^/^ J dk i2cuiy^^^ a (k)+ a t(-k) exp(ikx), (2) 

TT (x) = -t (27r)-'/' j dk K/2)'/' [a (k)- a t(-k)] exp(zkx), (3) 

^P (x) = i2n)-'/' I dp [m,/Eiy'' Yl U^ (P, A (p, r) exp(zpx). (4) 
In the last expression we have introduced the following matrix notations 





f (P.r) = '.'"•■ = o"""' . (6) 



In other words, the lower index i (the energy-sign index) in Ui (p, r)and (p, r)reflects the 
particle-antiparticle degrees of freedom in the Dirac formalism. The quantities k, p and r 

O O 

are the particle momenta and the fermion polarization index. The two Dirac spinors u and v 
satisfy the conventional equations — mo) u (p, r) = and + mo) v (p, r) = with 
^° = -E'p7° — P7 . The relativistic bare energies are expressed as = p^ + ml and uj^ = 

^/k^ + /ip, where the unknown values tuq and /io play role of the bare (nonrenormalized) 
masses. As in Ref. we use throughout this paper the definition of the Dirac matrices, 

with orthonormalization as in Bjorken-Drell (Ref. j39|) for u and v, 



o o E° 

Ui ^ (p, r) Uj (p, r') = -^6i.j6r,r'. (7) 

mo 

o o , ° of o ot 

The operators a (k) and a '(k), b {p,r) and b {p,r), d {p,r) and d {p,r) satisfy the 
usual commutation relations 
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5(k - k') 



ot 



d{p,r),d {p',r')} =6rr'S{p-p') 



a (k) , a (k'] 
b ip,r)X (p',r') 



which, in particular, yield 

Fi (p, r) , Fj ^ (p', r') | = 6^J5r,r'6{p - p'). 
In the BPR the corresponding unperturbed Hamiltonian is 



Ho{a) = I dkuj^ a ^ (k) a (k) + / dpE° 



of o ot O 

b (p,r) 6 (p,r) + d (p,r)rf(p,r) 



(8) 



(9) 



(10) 



By definition, the one-bare-fermion and one-bare-meson states \p,r)° and |k)° are eigen- 
states of Hq with eigenvalues E^ = ^/p'^ + ml and = A/k^ + /ig. These states are 

o o o o , o 

built of the bare vacuum Qq, |p,r)° =b '(p,t) Qq and |k)° =a T(k) Qq, respectively. 

O Q o o 

As usual, the no-bare-particle state f2o is destroyed by the operators a, b and d, i- e., 

o o o o o o 

a (k) f2o=6 (p, r) ^lo=d (p, r) fio= for Vk, p, r. 



B. The mass-changing BogoHubov-type unitary transformation 

Now, let us consider a set a = (a, a''', ...) of destruction and creation operators for particles 
with generic masses m and yU. If m and fi assume the physical values, the representation 
refers to "bare particles with physical masses" (see Ref. |l5|). By definition, the operators 
a enter in the following expansions of the same fields: 

(x) = (27r)"^/^ j dk (2LUk)"^^^ [a(k) + a\-k)] exp(ikx), (11) 
TT (x) = -i (27r)"^/^ J dk (cUk/2)^/^ [a(k) - a^(-k)] exp(ikx), (12) 
^ (x) = {2n)-'/' j dp {m/E^f (p, r) F, (p, r) exp(2px), (13) 

r,i 

where u (p, r) and v (p, r) are the Dirac spinors, which satisfy the equations 

— m)u (p, r) = 
(/» + m)t> (p, r) = 

with = -Ep7° — P7, -Ep = ^Jp^ + rn^ and c^k = i/k^ + yU^. The spinor column [/ (p, r) is 

o o 

composed of the spinors u (p, r) and f (— p,r) just as the column U (p, r) of u (p,r) and 
^ (— P,t). Similarly, the operator column F(p,r) is composed of the operators 6(p, r) and 

o o ot 

(P) just as the column F (p, r) of b (p, r) and (p, r). 
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Comparing the expressions ()lH ) - p3p and Q-©, we find the hnks 

a(k)+at(-k) a(k) + at(-k) 



-, Vk, 



(14) 
(15) 



(k)- a t(-k)) = (a(k) - at(-k)) Vk, 
E U^ (P, r) A (p, r) = J^Y.U. (P' ^) (P' ^) ' ^P- (16) 

V P r,2 VP r^i 

o 

Moreover, operators a are assumed to meet the same commutation rules as operators a 
do. Thus, the relations © are conserved as well. The orthogonality of U (p) follows from 

Eq. ((7j), after the substitution mo m {E° Ep). Since the operators a and a meet the 
same commutation relations, it is reasonable to look for a similarity (unitary) transformation 

a (k) = Ta (k) Tt, b (p, r) = Tb (p, r) Tt, (p, r) = Trf (p, r) Tt, (17) 

that connects them. We confine ourselves to the form T = Tmes ® Tferm, where the UT's 
T^es and Tferm ^ct in mesonic and fermionic sectors respectively and generate the linear 
relations: 

a (k) = Tmesa (k) T^,, = ci(k)a(k) + C2{k)a\-k), (18) 
with real functions Ci(k) and C2(k) and 



(p, r) = T;e.„^Fi (p, r) T)^,^ = ^ 0^,,;,- ,.(p)F,- (p, r') 



(19) 



where for a given p the c-number coefficients Oi^r-jyip) form a unitary 4x4 matrix 0(p). 

Evidently, the constraint — c| = 1 is necessary to ensure the commutation 
[a(k),at(k')] =5(k-k'). 

The condition (fTHj) can be fulfilled if one introduces the following ansatz j47| 



T 



exp 



^ y dkxfc (at (k) at (-k) - a (k) a (-k)) 



(20) 



where Xk is some real function of = |k| 
In fact, then we get 



a (k) = coshxfca(k) + sinhxfcat(— k). 



and Eq. dH)) holds if 



exp Xk 




(21) 
(22) 



I.e., 



cosh Xk 
sinh Xk 




+ 
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It can be readily verified that with the choice ()22p the condition p5|) is automatically 
satisfied. 

In regards to the fermionic sector, starting from Eq. (fT^ and omitting the discrete 
indices, we find 

^ „ ^^^^ 



'^f/(p)0(p) = ^-t/(p). 



Further, the orthogonality for the spinors JJ enables us to write 
0(P) = 



mrriQ ° ^ 



mnio I M T (p, r) M (p, r') {p,r)v {—p,r') 

^p^p \v ^ (— p, r) u (p, r') V ^ (— p, r) v (— p, r') 

The diagonal elements of the matrix (j2H) turn out to be equal to each other: 

u ^ (p, r) u (p, r') =v ^ (— p, r) f (— p, r') = 



(24) 



'(Ep + m)(E° +mo) 



mrriQ 



1 + 



(Ep + m)(E; + mo) 
while the off-diagonal elements differ from each other only by sign: 



u (p, r) V (— p, r') 



° t 



2' 



(^p + m)(E° +mo) 
mrriQ 



-p,rjM (p,r 
1 



E° + mo 



(0, r) p7t> (0, r') 



where (0, r) and (0, r') are the Dirac spinors in the rest frame and the usual Dirac matrix 
a 
-a 



7 is expressed as 



This allows to rewrite the unitary transformation for fermions in the following way 

Oip) = Aip)I + Bip)pT, 

where, along with the unit matrix /, we have introduced the 4x4 matrix 

"0 a' 
-a 

in the space of energy-sign and polarization indices. 
The functions A{p) and B{p) are equal to: 



(25) 



A{p) 



1 l{ Ep + m){E°+mo) 
EpEl 



1 + 



(Ep + m)(E° +mo) 



B{p) = -. 



l{Ep + m){E°+mo) 
EpEl 



E° + mo Ep 



m 
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and unitarity (0^0 = 00^ = /) leads to the relation 

(p) + p^B^ b) = 1. (26) 
Furthermore, as in the meson case, one could start from the ansatz 



Tferm — exp 



rfpFt (p) h (p) F (p) 



(27) 



where h (p) is an antihermitian 4x4 matrix. With such a form we get 

O (p) = exp [h (p)] , (28) 
assuming that h (p) = G (p) T and taking into account that 

exp [G (p) r] = cos (||G (p)||) + sin (||G (p)||) . 

Here ||G (p)|| denotes the length of the vector G (p). 

Comparing this with Eq. (j2SI), we can write G (p) = pQ {p) where the real function Q (p) 
is determined by 

pB (p) = sin [pQ [p)] . (29) 

After we have specified the transformation in question, the free Hamiltonian can be 
expressed in terms of the operators a, 

(30) 

where 

Hpia) = j dkcuka^ (k) a (k) + j dpE^ [b^ (p, r) b (p, r) + (p, r) d (p, r)] , (31) 

Mren,mes («) = ^^^^ / — [a^ (k) a (k) + a (k) a (-k) + H.c] , (32) 
Mren, ferm («) = rn{mo ' '^n) j ^ ^ F} (p, t') Mij (p, r'; p, r) Fj (p, r) , (33) 

^ r,r' 

with Mij (p, r'; p, r) = t/j (p, r') Uj (p, r), where Ui (p, r') (i = 1, 2) is the element of the row 
[u{p,r) ,v {—p,r)]. The use of the subscript "ren" associated with the term "renormaliza- 
tion" may be justified a posteriori (see Subsect. 2.3). 

The linear expressions ()18|1 and perform the transition from the initial particle rep- 
resentation with primary (bare) masses to the auxiliary particle representation with other 

o 

{e.g., physical) masses. In the latter we have a new vacuum state Qq = T* VLq that is 

o o o 

destroyed by the new operators a, 6, d just like the operators a, 6, d do the same with the 

o 

primary vacuum f2o- It is noteworthy to observe that the r.h.s. of Eq. ()30|) involves the 
terms that do not conserve the number of bare particles with new masses. Simultaneously, 
we have constructed one of possible unitarily equivalent Fock representations for the system 
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of noninteracting bosons (mesons) and fermions (nucleons and antinucleons) . In order to 

O 

differentiate this representation from tlie primary representation ot we will refer to it as the 
a representation. Its vacuum state IIq , the one-particle states a'^f^o, &^fio and (i^f^o together 
with other many-meson and many-fermion Hp eigenstates form the basis of the Fock space, 

spanned by the primary basis f2o; f^O; h f2o) d ^ Vt^-, etc. 

At this point, one should note that the transition ()18|) has much in common with the so- 
called cosh-sinh transformation for pair of boson operators, introduced by Bogoliubov ji^l 
in the theory of the ^He superfluidity. However, such a resemblance is rather mathematical 
since the Bogoliubov transformation and the UT Tmes have different scopes. As a matter 
of facts, the Bogoliubov transformation is related to the Hamiltonian of weakly interacting 
bosons to reduce it to diagonal form in the representation of quasiparticles (for a clear and 
simple discussion on this subject we refer to Chapter 30 of Ref. [41]). In that respect, the 
argument Xk is determined there with the help of a completely different physical condition. 
Instead, in the case of the transformation ()18|) we deal with free bosons and move to the 
opposite direction: from the diagonal form of Eq. (fTUj) . to the form (jHUj) in which the number 
of bosons with new mass fi is not conserved. 

Similarly, the transition (fT^ with the matrix O determined by (^3)) is an analog of the 
Bogoliubov u-v transformation in the theory of superconductivi ty ( see, e.g., Ref. |4^). The 
latter is a particular case of the general canonical transformation |4q for operators of fermion 
pairs, after Koppe and Miihlschlegel [43^ . 



C. Unitary clothing transformations and mass renormalization 

We discuss now some ke y p oints of the clothing procedure in QFT (more details can 
be found in Refs. and (s^). To be specific, we consider the model of meson- nucleon 

pseudo scalar (PS) coupling, in which Vq = Vo{a) = igQipj^tpcj) = ig^ F(a;mo,yUo)- We 
refer to models of this type as Yukawa models. At this point, one should stress that the 
described transition to the a-representation does not change the primary interaction Vq = 
igoY^a] m, n) with F(a; m, /x) = Y{a] mo, /xq)- 

By using the decomposition (jHUj) . we divide the total Hamiltonian H{a) into the new free 
part Hf^o) and the new interaction term Hj{a), 



H = H{a) = HF{a) + Hi{a), 

Hj{a) = V{a) + Mren,mes{a) + M,e„,/erm(a) + K-en(a), (34) 

where the operators Mren{c() are considered as the mass counterterms. Here we assume that 
the operator V = igtp'j^ilxf) depends on the "physical" value of the strength parameter g, so 
that a vertex counterterm appears in the interaction term Hj{a), Vren{,Oi) = Vo(a) — V{a). 
One could proceed further in the discussion by introducing properly regularized interaction 
V = Vreg with the cutoff vertex functions, but we will postpone such an extension to another 
occasion. 

In order to motivate our transition to the CPR, we recall again that the vacuum state Qq, 
with no bare particles with physical masses, and the one-particle states a"''(k)fio5 ^^(P5''")^o 
and d^{p,r)Qo are Hp eigenvectors, but are not eigenstates of the total Hamiltonian H. 
This is clearly seen in the current Yukawa model, where the interaction operator in the 
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a-representation looks like 

Vo{a) = J dkV^a{k) + H.c, 

V"" = J dp'dp F} (p', r') V^^ (p', r'; p, r) F, (p, r) . (35) 

r,r' 

In the above equation, we have introduced the 2x2 c-number matrices (cfr. App. A of 
Ref. il5,J), 

(p', r'- p, r) = -^,^^^—5 (P + k - p') (p', r') ^,U, (p, r) . (36) 

This interaction contains, for instance, terms responsible for the "non-diagonal" tran- 
sitions, such as VIq — > a^Wd^Vtf) and Ii^VIq — > a^b^bb^Q,o. Together with the a^a^ and 
b'^d'^ contributions to Mren,mes and Mrenjerm they prevent the aforementioned vectors to be 
eigenvectors of the total Hamiltonian H. 

To overcome this problem, the CPR introduces another representation of the total field 
Hamiltonian, 

H = K{a,) = Kpia,) + Ki{a,), (37) 

where the decomposition into a free part Kpiotc^ and an interaction part Ki{ac) depend on 
newly defined destruction and creation operators etc, 

a, (k) (at (k)) , 6, (p, r) (fot (p, r)) , 4 (p, r) (4 (p, r)) , Vk, p, r. (38) 

These are called the clothed particle operators|4^. By definition, they have the following 
properties: 

i) The physical vacuum (the H lowest eigenstate) must coincide with a new no-particle 
state f2, i.e., the state that obeys the equations 

a, (k) P) = b, (p, r) \n) = 4 (p, r) \n) = 0, Vk, p, r (39) 

= 1. 

ii) New one-clothed-particle states |k)c = a[ (k) Vt etc. are the eigenvectors of both Kp 
and K. 

iii) The spectrum of indices that enumerate the new operators must be the same as that 
for the bare ones. 

iv) The new operators satisfy the same commutation rules as do their bare counterparts. 
For instance, 

K (k) , at (k')] = 5(k - k'), [b, (p, r) , fot (p', r') ] = {4 (p, r) , dl (p', r') } = 5„,5(p - p'). 

(40) 

To be more specific, the property ii) implies that if the clothed meson state |k)c belongs 
to the eigenvalue tUk = ^/k^ + /^^ of the operator Kp, then K\k.)c = uj]^\k.)c. In other words, 
the operator Kj to be found must possess the property: 

Kj\k), = 0. (41) 
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The same is valid for the clothed fermion states. This property defines an important dis- 
tinctive feature of the CPR. 

Now, when finding the operators etc as functions of a, in order to meet the property iii) we 
suppose ae = W^aW, where W = W{a) = Wia^) = expR{ac) is a UT {WW^ = W^W = 
1). Before constructing its generator R{ac) = —R^{ac), let us rewrite the total Hamiltonian 
as 



H = H{a) = H{W (aj oicW^ {oic)) = W (a^) H (aj {a^) = K (a, 
= Hp («c) + Hi (ae) + [i?, Hp] + [i?, Hi] 

+ ^[i?, [i?, Hp]] + ^ [R, [R, Hi]] + ... (42) 

It is important to realize that the operator K {a^ is the same Hamiltonian as H{a) but 
it has another dependence on its argument etc as compared to H{a) because it refers to 
a different representation. Also one should note that the new free part Kp = Kp{ac) 7^ 
Hp{a) = Hp, but Kp{ac) = Hp{ac)- Hence, comparing Eqs. (jHTj) and we see that 
Ki (ttc) = K (aj - Hp (aj- 

Eq. ()42|) gives a practical recipe for the K (etc) calculation: at the beginning one replaces 
a by Oc in the initial expression H (a) and then calculates W (oc) H (oc) (oc) using Eqs. 
(j321) and (jlOI). The above transition H (a) — >• H (uc) generates a new operator H (ac) as 
compared to H (a), but Eq. (ji^ show that W (oc) H (a^) (etc) turns out to be equal to 
the original total Hamiltonian. 

Further, to meet the requirements i) and ii), the r.h.s. of Eq. (j42j) must not contain some 
undesirable terms that prevent the no-clothed-particle state fl and one-clothed-particle 
states to be eigenstates of the total Hamiltonian. Such terms (we call them bad as in Ref. 
[3^) enter in the operator V{ac) that is derived from V{a) by means of the replacement 
a — > Oc- By definition, these terms, taken together with their H.c. counterparts, do not 
destroy the physical vacuum and simultaneously the one-clothed-particle states. 

In case of the Yukawa model all terms in the r.h.s. of Eq. (j35|) are bad, and to eliminate 
them from K (etc) we choose an R that satisfies the relation 

V + [R, Hp]=0 . (43) 

It turns out that a solution of this equation can be represented as [l5l | 

POO 

R=-i lim / V{t)e-''dt , (44) 







where V(t) = exp(iHpt)Vexp{—iHpt) is the interaction operator in the Dirac picture. 
According to [15], the corresponding operator R can be expressed as R = TZ — TZ^ with 



4^ 



7^ = y c/k^^ac(k), 
j dp'dp Fl (p', r') Rl^ (p', r'; p, r) F,,, (p, r) . (45) 
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Here, unlike the fermion operators F and F'^ in Eq. ()35|) the operator column F^ and row 
Fl are composed of the clothed nucleon and antinucleon operators. As shown in f]^, the 
c-number 2x2 matrix i?^ is determined by 

Rl, (P'Z; p ,r) = V^^ (p'Z; p ,r) / [(-1)^"^ E^. - (-l)^'"^ E^ - u^] , (2, j = 1, 2) . (46) 

One should note that the solution of Eq. (j43|) exists if /i < 2m. This condition has a 
clear physical meaning, in that, if /i > 2m, the meson can decay into the //-pair making 
one-meson states not stable, namely, they cannot be H eigenvectors. 

Once [i?. Hp] = —V, Eq. can be rewritten as follows 

K (a,) = Hp (a,) + M^en ("c) + Ken M + ^ [i?, + [i?, M,e„] + ^ [i?, [i?, \/]] + ... . (47) 

Doing so, we have removed from K (etc) all the bad terms of the gf^-order and left the 
contributions up to the (7^-order terms assuming that Mren ~ Oi^g"^) and Vren ~ 0{g^). 

However, the r.h.s. of Eq. ()47p embodies other bad terms of the g"^- and higher orders. 
The commutator [i?, V] involves the terms g%'^^d\a\a\ , which do not destroy the physical 
vacuum VL . In addition, we find in [R, V] the terms g'^bld^alac, which neither destroy alfl, 
nor retain it with a multiplicative factor. These and similar bad terms can be eliminated 
with one more transformation, in a way analogous to the described above. 

The commutator [R, V] also contains the meson and fermion two-operator terms a^a, aa, 
b^b, b^d^ and so on, whose structure repeats the structure of the expansions and (jHS))- Not 
all of them are bad (for instance, a^a). It is required that the "diagonal" (particle-conserving 
number) species of the a^a-, b^b-, and d^d- types cancel the corresponding contributions to 
the mass counterterms Mren, mes («c) and Mren, ferm («c) • Note that it is sufficient to evaluate 
the mass shifts m — mo and fi^ — /Iq in the (yf^-order; since the same operator structure will 
appear in higher orders in the coupling constant we can extend this requirement to determine 
these mass shifts order by order. 

The first results in this direction have been obtained in for interacting pions 

and nucleons in the Yukawa model with the PS coupling. For example, according to |il5il] . 
the meson mass shift in the g"^ -order is equal to 

5fi' ^fil-fi' = ^ f^ll + ^ ] ■ (48) 



The nucleon mass shift evaluated in j36| in the same order can be written as 

g^ 

Mq / ni^ -pq + pq \ 

J E^X 2 [m2 -pq]-^2+ 2 [m^ + pq] - fi"" } ' 

Here we have adopted the four-vector notation, namely p = {Ep,p), k = (cjk, k) and q = 
(i^q,q)- 

We observe that, being expressed through the explicitly covariant integrals in the r.h.s 
of Eqs. PH| and these quantities do not depend on the particle momenta. It turns 
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out [36] that these integrals coincide with the correspond ing one-loop Feynman integrals. 
In this connection, we observe that recent calculations in |4^ have been carried out using 
Okubo's idea with a toy boson-fermion interaction model. The result obtained there for the 
fermion mass shift reproduces the corresponding expression deduced via Feynman-diagram 
technique. 

Up to now, the way the mass parameters have been introduced appears, to a great extent, 
quite artificial. However, once m and /i are identified with the observed physical masses, 
then, the eigenvalue Ep in the equation H |p, r) = Ep\p,r), with Ep = y/p^ + m^, becomes 
the fermion energy. And analogously, one can discuss the appearance in the CPR of the 
meson physical mass fi, which is different from the trial mass fiQ. This represents, in our 
opinion, a very natural way to introduce (renormalized) masses in the system. 

Finally, we would like to stress that the mass renormalization is implemented here in 
concomitance with the construction of relativistic interactions involved in operator Kj and 
describing the physical exchange-type processes between clothed particles. In the following, 
we will give a reasonable approximation to the operator Kj that must meet the condition 



III. THE CONSTRUCTION OF RELATIVISTIC INTERACTIONS VIA UNI- 
TARY CLOTHING TRANSFORMATIONS 

Before deriving analytic expressions for the interactions in the CPR one should emphasize 
that the corresponding operators stem primarily from the commutators [-R, l^], [R, [R,V]], 
etc. in the r.h.s. of Eq. (j47p . It is convenient to classify such operators by the numbers of 
creation and annihilation operators that appear in the normal ordered product. According 
to this classification an operator is of class [m,n\ if its normal product is made of 
m creation and n annihilation operators. In fact, after normal ordering of the creation 
(destruction) operators ac one can separate out the "good" four-operator contributions of 
the ^f^-order V^^l^ (ac) ~ a\h\achc., h\h\hchc and other operators of the class [2, 2] from the 
term | [i?, V], the "good" five-operator contributions of the -order V^H^ (ac) ~ a\h\h\hchc., 
a\h\h\hchc and other operators of the class [3,2] from the term | [i?, [i?, V^]], etc. Some of 
them will be interpreted below. 

In general, it is important to keep in mind that the so-called "good" terms and their H.c. 
counterparts, by definition, destroy the no-clothed-particle and one-clothed-particle states 
(cf. Eq. (jHJ). Unlike this, as mentioned in Subsect. 2.3, the bad terms have not such 
property, viz., even if one of them destroys such clothed states, its H.c. counterpart does 
not. Of course, among the operators involved in the total Hamiltonian H there are the 
operators of the class [1, 1] (for instance, ajoc) that ensure these states to be H eigenstates. 

As the next step of the clothing procedure we consider a recursive continuation of Eq. 
(gZj) with 

K («c) = Hp (aj + Vc (ttc) = W'K {a',) W'^ = K' {a',) = Hp (a^) + («c) , (50) 

where the operators a'^ in (jSUj) are related to the operators Uc in (jTTj) via the similar- 
ity transformation a'^ = W'^ acW . The interaction part Vc (etc) = H — Hp (oc) = 
W (etc) H (ttc) W"^ (oc) — Hp (etc) resulting from the first clothing can be represented as 



rc(«c) = exp(^)i^(i?; A) 



- Hp {ac 

A=0 
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where we follow the notation A{B; A) = e^^Ae~'^^ with the real parameter A for any oper- 
ators A and B. Obviously, A = A{B;0). According to this notation, the condition imposed 
on the generator R by Eq. (jlH)) reads 



0, 



A=0 



and the interaction generated by the first clothing UT acquires the expression 

d 



Vc = exp( — ) [MreniR; A) + VreniR] A)] 



+ 



A=0 



du 



exp(^— ) — exp(M— ) 



'dA d\ 
This compact expression can be expanded in series 

1 c/" 



V{R- A) 



(51) 



A=0 



Vc = Mren + Ken + [M,e„(i?; A) + VreniR; A) + /„\/(i?; A)] 



n=l 



A=0 



with the coefficients /„ = Keeping in mind the relationship 



[Ar 



d' 



dX^ 



[B, [B,...[B,A]... 



(52) 



A=0 



with n brackets (n = 1, 2, ...) it is obvious that ()5H) is a compact representation of Vc as the 
expansion in the multiple commutators [MrenT, [VrenT ^^nd [V]"'. This expansion does not 
contain any terms of the first order in the coupling constant. 

The second clothing UT W = exp [i?*^^)] has the generator R^'^^ which meets the condition 



v,':i+[R^^H,]=o 



or 



V^,+ ^H^{R^'^;X) 



d 



dX 



(53) 



A=0 



to remove from Vc all the g'^-order bad terms. The latter belong to the classes [2, 0], [2, 1], 
[3, 0], [3, 1] and [4, 0]. Their sum V^^l enters in the decomposition V^^^ = | [R, V] + Mien = 

(2) (2) (2) 

^good + Kad- "'■^ tum, V(j deuotes the lowest order contributions to the series in g, 



Vc = Vaood + V 



^bad 



n=2 



with 

Vgood = V^2d + ^iold + ^iood + terms of the ^(^-and higher orders 

and an analogous splitting of the operator Vc{bad), whose terms do not allow the no-clothed- 
particle and one-clothed-particle states to be H eigenstates. 
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Further, repeating the tricks that result in Eq. ()51|) for the after-first-clothing interactions 
and using the condition ()53|). we arrive to 

(a2 =H-Hf (a;) = K,„rf(fl'''; A) + Vi^_„,t+ 



1 rf" 



^,ood(i?('^ A) + Had,.e.*(i?('^ A) + /„V;S(/?('); A) 



(54) 



A=0 



where Had.rest = Vbad - V^feJ = 0{g^). Along with the framework interactions yjfj^ 
2,3,...) expressed through the new clothed operators a'^ this equation enables us to evaluate 
the ^f'^-order corrections (n = 4, 5, . . . ) to them. As an illustration, omitting the contribu- 
tions of the - and higher orders we find 



(n) 



(3) 



good 



+ k2 + 



good 



- good 2 - 



' ^bad 



~^ ^good 



a + 



(55) 



The r.h.s. of this equation involves the good interaction operators of primary interest since 
they can be associated with different processes in the meson-nucleon systems. In particular, 

(2) 

as aforementioned, the term V^^^^ consists of the interactions of the class [2,2]. In order to 
find the corrections of the (/^-order to them one needs to separate out the operators of the 
same class of the three last terms in the r.h.s. of We have pointed out the g^-oider 

(3) 

contribution V^^j as the first candidate for eliminating via the third clothing UT. 



A. Details of Calculations 

We have seen how the interactions between clothed (physical) particles can be constructed 
when handling the multiple commutators [V]"' (n=l,2,...) with n-brackets. In the framework 
of the Yukawa model or any other field model with a polynomial interaction the operator V 
can be represented in the following symbolic form: 

V = f *m + H.c, (56) 

where / * m is a polynomial composed of products of fermionic and mesonic operators. 

In general, to obtain recursive relations for the commutators of increasing complexity, it 
is convenient to write in accordance with Eq.()52p 

[Vr = lim 4^ (e^^/ * me^^^) + H.c. . (57) 

Taking into account that R is an antihermitian operator, we have: 

[Vr = hm^{f{X)*m{X)) + H.c.. (58) 

For brevity, we assume /(-R; A) = /(A) and m{R; A) = m(A). 
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Then, using the Leibnitz formula 



we find 



whence 



dX 



- (/(A) * m(A)) = C'nifWr' * [miX)Y 



s=0 



lim ^ (/(A) * m(A)) = J] Q[/(0)]— * [m{0)Y 



s=0 



s=0 



At the first stage of our procedure for the Yukawa model 



(59) 



Vfe„^ = V (aj = / dk/(k)m(k) + H.c, R = Ri (aj = / dkR^ m(k) - H.c, (60) 



with /(k) = and m(k) = ac(k). For such V^ad formula where [/]" * * [m]* = 

/ c/k[/(k)]"-^[m(k)]^ yields 

[Had]' = l^ik [[/(k)]^m(k) + /(k)[m(k)]i] + H.c, (61) 

[Had]' = y"rfk [[/(k)]2m(k) + 2[/(k)]i[m(k)]i + /(k)[m(k)]2] + H.c, (62) 
From these equations it follows that 



[R^V] = [VY = / dkidka 



T>ki 



a(k2)a(ki 



^(k2)a(ki) + Hi^ii?''^t5(k^ - ka)} + i^.c. 



(63) 



[i?, [i?, V]] = [Vf = j dk.dk^dks {A^ik^, ks, k3)at(k2)at(ki)a(k3) 

+A2{ki, ka, k3)a^(k2)a(ki)a(k3) + A3(ki, k^, k3)a(k2)a(ki)a(k3) 
+A,{k,, k2)at(k2)5(ki - ks) + A5(ki, k2)a(k2)5(ki - k3) } + H.c, 



Al(ki,k2,k3) = R^\ 
v43(ki,k2,k3) = 



^ki yk2 



,yl2(ki,k2,k3) 



R^\ 



R^' 



^ki yk2 



A,ik,,k,) = R^^ 



R^\V^^^ 



+ 



t 



+ 



^k2 p.ki 
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A5(ki,k2) = B^' 

+2R 



+ 



V 



kit 



ki 



^k2 T^kit 



+ 



^k2 J>ki 



^kit + yki 



(64) 



For brevity, the superscript c has been omitted in the r.h.s. of Eqs. (jUHj) and 

Performing the normal ordering of the fermion operators in Eqs. and (jMj) . we get a 
simple recipe to select the 2 < — 2 and 2 < — ^ 3 interaction operators of the g^- and (yf^-orders 
between the partially clothed pions, nucleon and antinucleons (in particular, vrA^ tcN, 
NN —>■ NN and NN irNN). At the same time this algebraic technique enables to 
select the two-operator (one-body) contributions which cancel the meson and fermion mass 
counterterms Mren{mes) and Mreniferm) in the ^f^-order (details are in [l^ and [36]). In ad- 
dition, we encounter the three-operator (vertex-like) "radiative" corrections which together 
with the similar terms from the commutators [R, Mreni'm^es)] and [R, Mren{fGf"m)] cancel 
the "charge" counterterm Vren in the ^f'^-order. The remaining bad terms must be removed 
via successive clothing UT's. 

Along this guideline, we arrive at the decomposition: 



= Kpiac) + K{NN NN) + K{NN NN) + K{NN NN) 
+ K{ttN ttN) + K{nN nN) + K{nn ^ NN) 
+ K{NN ^ ttNN) + K{NN ^ nNN) + K{NN ^ nNN) 
+ K{NN ^ tttttt) + K{ttN ^ ttttN) + K{-nN ^ vrvriV) + . . . 



(65) 



where the interactions between the clothed nucleons (A^), antinucleons (A^) and pions (tt) 
have been separated out. 



B. Pion-nucleon interaction operator 

To obtain the explicit expression for the operator K{7cN — > vrA^) one needs to separate 



out the blbc-tjpe terms form the commutator 



^k2t T>ki 



and its H.c. part in the r.h.s. of 



Eq. (|63|) . The result can be represented in the following covariant (Feynman-like) form: 



A'(7rA^ vrA^) = j dpidparfkirfka Kjv(k2, P2; ki, pi)a1;(k2)6j(p2)ac(ki)6c(Pi) 



KAr(k2,P2;ki,pi) 



m 



2(27r)3 ^^k,^k2^p,^p2 



--S{pi + ki - p2 - k2) 



XU{P2) 



+ 



+ 



+ 



Henceforth the spin indices are omitted and all summations over them are implied. 



m(Pi). 

(66) 
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The ttN quasipotential being the kernel of an approximate integral equation (cf. Eq. 
4.9) in for the tiN scattering wave function in momentum space is determined as in 
3: 



Kiv(k2,P2;ki,pi) = (at(k2)6t(p2)n|ir(vriV ^ 7riV)|at(ki)6t(Pi)^^> 



(67) 



and it turns out to be equal to K-Ar(k2, p2; ki, pi). 

In order to interpret Eq. . we write an intermediate analytical result that leads to it, 
namely 



KAr(k2, P2; ki, Pi) 



m 



=5(pi + ki - p2 - k2) 



P2 



X 2"(P2)75 {P (k2, P2; ki, pi) + P (ki, pi; k2, P2)} 75m(pi) 



(68) 



P(k2,p2;ki,pi) = — 



A+ (g) 



+ 



A^- (g) 



A+ {q') 



AL (g') 



P2 



-Eq' — u;ki -Epj + -Eq/ — u;ki 



(69) 



with the standard notations A_|_ (A_) for the projection operators on the fermion positive 
(negative)-energy states A-t (g) = (± ^ + m) /2m. 

Separate contributions to the r.h.s. of Eq. (jU^ can be represented as in Fig.l, viz., via 
graphs a and 6, which correspond to the two terms in the first square brackets, and c and 
d associated with the two ones in the second square brackets. The translational invariance 
imposes the constraints ki + pi = q = k2 + P2 and pi — k2 = q' = P2 — ki upon the momenta 
involved. In other words, it implies that three-momentum is conserved in the each vertex of 
these graphs. 



(a) : k 




FIG. 1: Different contributions to the vrA^ quasipotential within the Old-Fashioned Perturbation 
Theory (OFPT). 
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Such graphs are well-known from the old-fashioned perturbation theory (OFPT), where, 
for example, the vrA^ scattering T-matrix in the ^f^-order is determined by the elements 

(k2,p2|T(2)(E + zO)m,pi) = {k2,P2\V{E + tO-Ho)''V\k„p,) 

= J2 (k2, P2\V\t) {E + tO- E,y' i^\VW, Pi) . 

i 

Here E the collision energy and the summation is performed over all permissible intermediate 
states (the Hq eigenstates) with energies Ei. In this context, the inverse energy denominators 
in the r.h.s. of Eq. ()69|) have the form {E—Ei)~^ with the appropriate values of E and E^. For 
example, (-Eq + E^^ + cukj"^ and {E^^ + Eq; — t^kj"^ may be related to the noncovariant 
propagators 

D'^ {E) \e=Ep^+w^,^ = {E - Ep^ - cjki - Ep.^ - cjk2 - Eq) ^ \E=Ep^+iOi,^ 
= -{Eq + Ep,+u;i,,)-\ 

and 

D (E) |£;=i?p^+a;kj = {E - Ep^ - Ep^ - Eq>) ^ \e=Ep^+uj^^ 
= - iEp2 + Eq' - t^ki) ^ , 

being associated with the five and three internal lines putting in between the dotted verticals, 
respectively, in graphs b and d. Both of them contain the NN pair in the intermediate states. 
Note that we have to distinguish the intermediate states of clothed particles in our approach 
from the intermediate states of bare particles with physical masses that usually appear in 
OFPT. 

In addition, we would like to note that the graphs in Fig. 1 are topologically equivalent to 
the four time-ordered Feynman diagrams displayed in Fig. 20 in Chapter 13 of Schweber's 
book [2^. However, in the Schrodinger picture used here, where all events are related to one 
and the same instant t = the use of such an analogy seems to be misleading. In fact, the 
line directions in Fig. 1 are given with the sole scope to discriminate between the nucleon 
and antinucleon states. 

The Feynman-like propagators in Eq. (jfifij) appear after trivial transformations, e.g., 
adding the contribution a and b, it is readily seen that 

1 

— m' 

where we have denoted the "left" Mandelstam vector si = {Ep^ + c^k^, pi + ki) = pi + ki. 
Now, taking into account the property 757^^75 = — 7/^, one shortly derives Eq. (jHEI)- 

We emphasize that energy conservation is not assumed in constructing this and other 
quasipotentials, and this forces us to give much more detailed and complicated expressions 
for these quasi-potentials. Indeed, the coefficients for the ajfejac&c-terms appearing in the 
K{ttN ttN) expansion generally do not fulfill the on-energy-shell condition 

-Epi + ^^ki = -Ep2 + t^ka, (70) 

In this connection, the "left" four- vector si is not necessarily equal to the "right" Mandel- 
stam vector S2 = P2 + k2- Clearly, if the condition (f7n|) is fulfilled, then pi — /c2 = u, so 



m 

e: 



A+(g) 



-^q ~ Ep^ — tUk^ 



+ 



Al (g) 

Eq + Ep^ + CJki 
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that the first (second) half-sum of the covariant propagators in the r.h.s. of Eq. ()66p is 
converted into the s-pole (w-pole) propagator typical of Feynman technique. Thus, only the 
corresponding on-energy-shell part of the vrA^ quasipotential is represented in Fig. 2 by the 
genuine Feynman diagram. 



'^1 



s-pole P2 P\ u-pole ^ 



FIG. 2: On-energy-shell contributions to the vriV-quasipotential. 

The quasipotential is nonlocal since the vertex factors and propagators in Eq. (jM)) are 
dependent not only on the relative three-momenta involved but also on their total three- 
momentum. Besides, this embodies the nonstatic (recoil) effects in all orders of the so-called 
expansion (see (45^). 

In the static limit, one can show that the off-energy- shell contributions differ from the on- 
energy-shell results by the values of the order (\1/rf where is the relative vriV momentum 
and Tj = rrifi/ {m + jj) is the reduced tiN mass. One needs additional investigations to show 
to what extent these differences can be neglected when describing vrA^ scattering . 

To our knowledge the formula (jHTjl (albeit given in other form) has been first derived in 
[3^ within the clothed particle approach. 

The quasipotential determined by Eq. has the same Feynman-like structure as the 
effective t^N interaction obtained in ji^ (see the final formula below Eq. (A. 8) therein). As 
mentioned in Introduction, the authors of ji^ have used the other realization of UT method, 
which is similar to the Frohlich transformation (see, e.g., 0|) in the theory of metals. 

In addition, one can show that the coefficients VJrAr(k2, P2; ki, pi) in Eq. (jMj) do not 
differ from those obtained in Ref. [13 (after some simplifications of the model Hamiltonian 
used there). However, the corresponding pion-nucleon interaction operator from Ref. 
contains the bare particle a'^fe'^afe-terms instead of the clothed particle a\h\achc-o\ies in our 
expression for K{ttN — > vrA^). It would be interesting and worthwhile to check if the same 
coincidence of results will be obtained for higher-order terms in the coupling constant. 



C. Nucleon-nucleon interaction operator 

After normal ordering of the fermion operators in the third integral in the r.h.s. of 
Eq. (j63j) and its Hermitian conjugate part the A^A^ — > A^A^ interaction operator can be 
represented in the form: 



K{NN^NN) = / rfpidp2rfpWp'2Viv^(p;,p'2;Pi,P2)feKp'i)&KP2)^c(Pi)&c(P2), 
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1 (j^ 1 TT? 
Viviv(p'l, P2; Pi, P2) = -1^7:^1^, 7====5(p'l + P2 - Pi - P2 



xm(p'i)75m(Pi) \ ^ + ^ \ m(P2)75m(p2) (71) 

As in case of the vrA^, we encounter here the noncovariant ( "nonrelativistic" ) propagators 
hnear in the intermediate energies. It allows us to relate the two graphs in Fig. 3 to the 
r.h.s. of Eq. 
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FIG. 3: Illustration of Eq. (f7T|) as a sum of two OFPT diagrams with the intermediate pion (dashed 
lines) on its mass shell. 

The left and right diagrams in Fig. 3 are associated with the OFPT denominators, 

D (E) \e=e^^+Ep^ = [E- Ep'^ - Ep^ - t^p^-pj \e=e^^+e^^ 



Ep, - Ep>^ - ^p;-pi 



and 



D (E) \e=e,+e, = [E- Ep^ - Ep'^ - ujp'^-p^) \e=e,+e. 



= Ep'^ - Ep, - iOp^-p'^, 

respectively (cf. discussion in Sect. 3.2). 
Though Eq.(f7T|) is equivalent to 

1 g'^ ttP 

Vnn{p'i, P2; Pi, P2) = -9 /9 /P F F F "^(Pl + P2 - Pi - P2) 

2 {27ry ^ Ep^Ep^Ep,^Ep,^ 

xm(p;)75m(pi)- ]— ^m(p2)75m(p2), (72) 

the occurrence of such primary denominators is typical of the three-dimensional formalism 
developed here. 

The corresponding relativistic and properly symmetrized A^A^ interaction, is given by the 
quasipotential 

V)v7v(p'i,P2;Pi,P2) = (6l(p;)6l(p;,)^^|ir(A^A^-.A^A^)|6t(p^)6t(p^)^^) 

= 2 {^JV7v(pi, P2; Pi, P2) - Vnn{v'i, P2; P2, Pi) 

-V^7V(P2, pi; Pi, P2) + 1^iViv(P2, pi; P2, Pi)} , (73) 
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(see formula (4.16) from Ref. |l5l|). 

Substituting (f7^ into (fTHj). we express the quasipotential of interest through the covariant 
(Feynman-hke) "propagators" , 

VMp[, Pi; Pi. P.) = - 2 p^5-^gpp=^*(p; + p;, - P, - p.) 

xs(p;)7.«(p0i + "(p'Mp^) - (1 " 2)- (74) 

The r.h.s. of this equation consists of the "direct" term written exphcitly and the "ex- 
change" one (1 2) with the prescription that the variables with label 1 and 2 should be 
mutually interchanged. Formula ()74|1 determines the A^A^ part of an one-boson-exchange 
interaction derived via the Okubo transformation method in |18j] (see also ji^l ) taking into 
account the pion and heavy-meson exchanges. 

As has been emphasized in Ref. [18], a distinctive feature of potential ()74|) is the presence 
of a covariant (Feynman-like) "propagator" 

1 r 1 1 1 

+ 7:— • (75) 



2 {{pi- p'lY - {p2 - ^2)^ - At^ J 
On the energy shell for the A^A^ scattering, that is 

Ei = Ep^ + Ep,^ = Epi^ + Epi^ = Ef, (76) 

this expression is converted into the genuine Feynman propagator which occurs when eval- 
uating the S - matrix for A^A^ scattering in the (7^-order. The corresponding graphs are 
displayed in Fig. 4. The quasipotential Vat at can be related to these Feynman diagrams, 
being different from the corresponding Feynman A^A^-scattering amplitude in two respects, 
namely, V^n does not contain the factor S{Ep-^ + Ep^ — Ep'^ — Ep'^) and "propagator" (f7^ 
is now related to the meson line in the left graph in Fig. 4. 
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FIG. 4: On-energy-shell contributions to the A^A^ interaction (the g^-ovdei Feynman diagrams for 
A^A^ scattering). 

One could take this analysis a little further, and show that the off-energy-shell deviations 
from the conventional on-energy-shell result obtained in the static approximation are of the 
order p^/m^, where is the relative A^A^ momentum. In other words, the corresponding 
off-energy-shell corrections have a relativistic origin. 

Explicit formulae for the A^A^-quasipotential generated by the exchange of heavier mesons 
are given in the next section. 
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D. A^A^ ttNN transition operators 



To separate the NN ^ nNN interaction operators one needs to perform the normal 
ordering of fermion operators in the fourth and fifth integrals of Eq. (jMj) and their Hermitian 
conjugate parts. After a lengthy algebra we get the pion production operator, 



K{NN nNN) 



J dpidp2dp[dp2dkV^NNip'i, p'2, k; Pi, P2) 



xal{k)bM)bMMPiMP2), (77) 



V^NN (p'l, P2, k; Pi, P2) = V^NN (Feynman - like) + V^nn (off - energy - shell) 



where 



K- AT AT (Feynman — like) 



m 



n(p^)75^(p2) , , 

{P2-P2) - 



1 1 



^(Pi + P2 - pi - P2 - k) 



'^1+ /t + m pi— }t + m 



u{Pi) 



(78) 



and 



K-Ar7v(off — energy — shell) = —i- ^ 



m 



U{P2)15U{P2] 



2uj, 



P2-P2 



3 (27r)9/2 ^2iu^Ep,Ep,Ep.Ep.^ 
1 1 



<5(Pl+P2-pi-P2-k) 



^P2 - Ep'^ +^P2'P'2 



xu{p[) 



+ 



1 



^1+ /c + m pi— + m 

1 



u[Pi] 



Ep2 - Ep'^ - ujp^-p'^ Ep>^ + c^k - Ep^ - Ljp.^-p'^ 
1 1 



+ 



1 



/c + m pi— + ra 
1 



^IPiJ 



P2 - ^P2-pi. -^P'l + ^^k - Ep^ - UJp.^-p'^ 



Ep^ — Er,l_ — CO- 

X^(Pi 



1 1 

+ 



u{Pi, 



1 



Eoo — E^' + 



.p2 ^p. , -'pa-p^ Ep'^ +UJk- -Epi + ^p2-p^ 



X ^^(Pl. 



+ 



{+j^- + m pi- /i- + m 



u[Pi, 



(79) 
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with q = P2 — P2 and g_ = (cjq, — q). 
We have used the relation 



[{x-yf-l/] ^ = {[Xo - VO - UJ^-y] ^ - [Xo - yo + ^K-y] ^) , (80) 



UJx-y 



that holds for any four-vectors x = (xo,x) and y = {yo,y) with lOz = ^/z"^ + 
On the energy shell: 

Ep, + Ep, = Ep. + Ep. + (81) 

expression ()78|) yields the pion production T matrix in the g^-order, i.e., to the first nonva- 
nishing approximation in the CPR. We remind that a clothed pion cannot be absorbed or 
emitted by a clothed nucleon. At the same time the r.h.s. of ()79|) becomes equal to zero 
since under condition (|HT] ) the propagators in its round brackets cancel each other. 

By analogy with the interaction between clothed nucleons we could introduce a quasipo- 
tential defined by the matrix elements of the operator (j77j) between the properly symmetrized 
clothed NN and vrA^A^ states 

KiV7v(p'i,P2,k;pi,P2) = (at(k)6t(p;)6t(p'^)(^|;^(iViV ^ ^rNN)\bl{p^)blip2)n) 

= - {Kjviv(pi, P2, k; pi, P2) - V^NNip'i, P2, k; p2, pi) 

-K7VJv(P2, pi, k; pi, P2) + K7VJv(P2, pi, k; P2, Pi)} (82) 

Of course, such matrix elements characterize only some part of that physical input which 
contributes to pion- nuclear dynamics with the operator ()77|) included in the interaction Kj. 

All the terms involved in the K-atat quasipotential may be divided into the two groups. 
The first group refers to the pion production mechanisms shown schematically in Fig. 5. 

Note that the three left graphs in Fig. 5 exemplify the processes in which the pion created 
by one nucleon is scattered by another one, while the three graphs on the r.h.s. show the 
processes, where one nucleon emits two pions one of which is absorbed by another nucleon. 
Again, on the energy shell the sum of all these contributions can be represented by the 
Feynman graph (Fig. 6, left) with the pion exchange preceding the pion emission. 

The second group of terms which could be extracted from Eqs. ()78j) and (f7n|l are shown 
in Fig. 7. 

On the energy shell, the sum of these contributions can be represented by the Feynman 
graph (Fig. 6, right) with the pion emission preceding the pion exchange. Therefore, in the 
first nonvanishing approximation a pion can be absorbed/emitted only by a correlated NN 
pair through a mechanism of the ^f^-order. 

In Ref. j2J] the UT method has been applied to the same A^A^ iiNN process assuming 
a meson-nucleon PV-coupling. Taking into account this difference in the coupling scheme 
and relying upon our Eqs. (f75|) and (fTUI). one can reproduce those effective interactions 
given first in 
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In search of physical interpretation of the obtained results let us consider the expression 

^(P2)75«(P2) 



- Ep'^ - UJp2~p'2 



1 1 

')[+ /c + m pi— + m 



w(pi). (83) 



We could write 

[Ep^ - Ep/ - a;p2_p/] ^ = [E-Ep^- Ep/ - ujp^-p'^ 
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FIG. 5: Illustration of the "retarded" pion production mechanisms on the A^A^ pair in the g^-order. 
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FIG. 6: The conventional Feynman diagrams for the A^A^ — > ttNN process in the (7^~order. 



with E = Ep^ + Ep2 that is prompted with the graphical representation of ()83p (more 
exactly, its part) in Fig. 8a, where all vertices displayed by circles ("interaction points") are 
connected by bold lines (among them the nucleon line has the right arrow to distinguish 
nucleons and antinucleons). 

As admitted in the OFPT, the particle three-momenta are related to the lines in Fig. 8a 
via the following prescription: when moving from the left to the right, a sum of "outgoing" 
momenta is equal to a sum of "ingoing" ones. Here the term "ingoing" ("outgoing") is 
referred to the corresponding line or lines lying on the left (right) from the dotted vertical 
(a "phantom" line) that passes through a given vertex. These vertices and phantoms are 
enumerated in the alphabetical order 1, 2, 3, ... Accordingly, we have 



P2 = e + P2, pi + e = q = k + pi. 



(84) 
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FIG. 7: Illustration of the "advanced" pion production mechanisms on the A'^A'^ pair in the g^- 
order. 
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FIG. 8: Graphical representation of Eqs. H83|) . 

Evidently, being taken together, these equations provide the total momentum conservation, 

pi + P2 = k + p'i + P2. (85) 

In parallel, we could consider Fig. 8b, where p2 = e+p2, e+pi+q = and q+p'^^ + k = 0. 
In particular, it means that — q = p'^ + k = e + pi = p2 — P2 in accordance with (jHH|) . 

As in case of Fig. 8a, the three internal lines between the phantoms 1 and 2 can be 
associated with the OFFT propagator, \Ep.-^ — E^i^ — cUpj.p/,] 

Regarding the lines between the phantoms 2 and 3 in Fig. 8a and Fig. 8b one can make 
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up the propagators (cf. discussion in Sect. 3.2) 



and 

- - Ep> - - Ep, - Ep, - ^p,_p. ] ' = - [Ep; +UJk + Eq] ' (87) 

with E = Ep^ + Ep'^ + ujp2-p'^ . 

Handhng the off-energy-shell graphs such as in Figs. 8a and 8b we cannot a priori ignore 
another option: with E = Ep^ +-^P2 +^p2-p'2 f|86|) and E = Ep'^ +-^P2 +'^k i^i (IHZ|)- Both 
options contribute to the r.h.s. of (jZHI), being equal to each other on the energy shell. 

Doing so and using the well-known rules, where, for instance, the vertices with the lable 
2 in Fig. 8a and Fig. 8b are equivalent to the factors 

=n(q)75«(pi) (88) 



2(27r)^a;p,_p, V^pi^q 



m 



and 

9 

a/2 (27r)a;k V^p'^E^ 
respectively, it is easy to get the matrix element 



^(P'i)75w(q), (89) 



u 



[p[)[^[+ /c + m]-\{q) (90) 



typical of the single- nucleon matrix elements involved in (|78p and (|79|) . 

Relying upon the discussion below (jTO)) one may say that the quantity is related 
to the "right-hand" s-pole mechanism of the off-energy-shell scattering of the intermediate 
meson (pion), that has been emitted by nucleon 2, on nucleon 1. 

As in Subsects. 3.2-3.3, one has to stress again that, unlike the Feynman diagrams in 
Fig. 6, the meson is on its mass shell. Being a mediator in forming different pion production 
mechanisms on the correlated pair of nucleons, it is intimately embedded in the diagram 
of Fig. 9a, where, first, one of the interacting nucleons "shakes off" the pion (the bright 
spot) which is absorbed, then, by the other nucleon (the hatched spot) with the subsequent 
("retarded") emission of the detected meson (pion). 




FIG. 9: Graphical representation for the "retarded" (a) and "advanced" (b) pion production 
mechanisms on the A^A^ pair. 

Of course, such "chronological" expressions which employ the terms "first" and "then" 
has nothing common with any real time development of the process NN ^ TfNN. In 
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fact, we imply a left-to-right-hand alignment of intermediate pion-fermion states that occur 
when evaluating the matrix elements (irNNlVln) {n\V\m) (m|V^|A^A^), where |m) and \n) are 
permissible Kp eigenstates. In this respect, it is useful to keep in mind the on-energy-shell 
relationship 

^ttA^A^I^ [R, [R, V]] INn"^ = {irNNlVGp (E) VGp (E) V\NN) 
with the resolvent expressed by 

Gf (E) = {E- KfY^ 

E = Ep^ + Ep2 = Epi^ + Epi^ + uj\^. 

Nevertheless, it is convenient to refer to the reaction mechanisms represented in Fig. 9a as 
"retarded", to differentiate them from the "advanced" ones displayed in Fig. 9b. 



IV. INCORPORATION OF HEAVY-MESON EXCHANGES 

In this section we have collected our more general evaluations of quasipotentials for the 
A^A^ — > A^A^ and A^A^ ttNN interactions that are due to the exchange of 6 = vr, rj, p, uj, S 
and a mesons. 



A. Pion and heavier-meson exchanges in NN NN process 

Here we consider the well-known Yukawa- type meson-nucleon couplings (see, e.g., [2ll |) 
additively involved in the primary total Hamiltonian. Application of our approach in this 
case results in the following expression for the A^A^ NN quasipotential 

~ 1 771^ 

Vnn{p'i, P2, Pi, P2) = - ^.^^^3 ^ „ „ „ ^(p'l + P2 - Pi - P2) 

2(27r)^ V^Pi^P2^pi^p'2 

X$^V^i^7v(Pi,P2;Pl,P2)- (91) 

b 

The contribution to this formula from the 6- and a- mesons exchange is given by 

V'' (W W r. r. ^ - n^h ^(p1)^(Pi)^(P2)^(P2) 
>/jVivlPl, P2, Pl> P2j - 9brib 7- -7T-2 -2 

\Pi Pi) n 

+ (1 ^ 2, 1' ^ 2') -(1^2)- (1' ^2'),b = 6, (7, hs = Ti- T2, K = 1. (92) 

As usually, the Pauli vector r is the matrix in the nucleon isospin space. The account for 
the TT- and rj- mesons exchange gives rise to the following contribution 

T/fc ^r.' r.'-r. r.^ ,2^ ^(P'l r'l hMPl ^iMp'^ r^)75W(P2 ^2) 
^Nn[Pij P2, Pi, P2j — 



+ (1^2,l'^2')-(1^2)-(l'^2'),6 = 7r,r7, = n ■ r2, h, = I- (93) 



30 



At last, the NN quasipotential originating from the p- and u- mesons exchange can be 
written as 



^iViv(p'i,P2;Pi,P2) = 9bhb 



u 



;pi)7M^(pi)Ar(pi - p'i)u{v'2)iMv>2) 



(Pl-p'l)2_^2 

+(1 ^ 2, 1' ^ 2') -(1^2)- (1' ^ 2'), 



(94) 



B. Heavy- meson exchanges in NN ^ irNN processes 

In regards to the NN ttNN process the corresponding quasipotential has the form 

im^ 5(pi + p2-p;-p'2-k) 



VnNN (Pi r[, P2 k; pi n, p2 r2) 



P2 



X X] ^=K?iviv(Pi P2 k; Pi ri, p2 r2). 

The (5- and cr- mesons exchange contribution to this formula is given by 
K?iViv(p'i'P2>k;pi,P2) = 9n9lu{v'2)u{p2)u{Pi) 



(95) 



X 



ti 



6,1 



''fe,2 



2a;'' 

q 



+ 



/; + m ^1— Jk — m 

''6,1 



£;p, - £;p. - uj\ E^^ - E^,^ + - c^i^ 

1 1 



+ 



''6,2 



''6,2 



''6,1 



X75^i(pi ri) + (1 ^ 2, 1' ^ 2') -(1^2)- (l' ^ 2'), 



The TT- and 77- mesons exchange gives the contribution 

K^7V7v(p'i ^1' P2 ^^2^ k; Pi ri, p2 r2) = g^g^uip'^ r2)75M(P2 r2)M(pi r'J 



(96) 



X 



+ 



1 1 

+ 



^6,1 



+ 



h,2 



ii+/; + m ^i— /c + m 



''6,1 



''6,2 



31 



1 



X'u(pi) + (1 ^ 2,1'^ 2') - (1 



fi 



2)-(l'-20, 



b^Tr,rj, = T1T2, ti^2 = T2TiTf , i^^i = i^_2 = 1^- (97) 

Finally, the ttNN quasipotential originating from the p- and uj- mesons exchange is 
determined by 

^i'NNip'i r'i,p'2 ^2, k; Pi ri, p2 ra) = 9n9lu{p2 r'2hMP2 r2)K''{p2 - P2)u{p'i r'l) 



X 



+ 



+ 



2^6,2(Pl-fe)^ 

/; + m ' pi—Jk-\-m {pi — — m? 
1 

+ 



^P2 ^P2 



X 



ti 



'h+ /i + m /i + m (p'l — qy — rn^ 



2uj^ 

q 



+ 



X 



^6,27^ 



<i7. 2tUp', + q-). 



^1+ ^_ + m j6i— + m {p'l + — 
X75^i(pi ri) + (1 ^ 2, 1' ^ 2') -(1^2)- (l' ^ 2'), 



iq'qyd)- (98) 



V. SUMMARY 



We have presented a field-theoretical method to construct in a consistent way relativistic 
(Hamiltonian) interactions for the meson-nucleon systems. In particular, the interaction op- 
erators for processes of the type NN — >■ NN, ttN — >■ ttN, and NN ttNN are derived on 
one and the same physical footing. The method is based on the unitary clothing approach 
which introduces a new representation of the primary total Hamiltonian in terms of the op- 
erators for creation and destruction of the so-called clothed particles, namely those particles 
that can be observed. Within the approach all interactions constructed are responsible for 
physical (not virtual) processes in a given system of interacting fields. Such interactions 
are Hermitian and energy independent including the off-energy-shell and recoil effects (the 
latter in all orders of the - expansion). The persistent clouds of virtual particles are 
no longer explicitly contained in the representation, and their effects are included in the 
properties of the clothed particles and in the interactions between them. We have also con- 
sidered the additional contributions to such interactions arising from heavy-meson exchange 
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mechanisms. Such mechanisms are an important subject of investigation in current research 
on few-body systems, and their true role in the pion-production processes (as well as when 
constructing three-nucleon forces ) has not yet been fully established. Such issues can be 
best investigated, in our opinion, via the unitary clothing transformation method which 
provides a well-defined, unambiguous definition and construction of relativistic interactions 
between particles with observable properties. 
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Such a form has been prompted by M. Shirokov during a tentative investigation carried out 
by him and one of us (A.S.) 

Here the term canonical applies to the transformation which does not violate the commutation 
relations for the fermions in question. 

As noted in Introduction, we employ the term "clothed" instead its synonymous "dressed", 
since the latter is sometimes used in a sense which differs from that defined below by the items 
i) - iv) 
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